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ABSTRACT 

Aims. We present the main results of an imaging survey of possible young massive clusters (YMC) in M31 performed with the Wide 
Field and Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST), with the aim of estimating their age and their mass. 
We obtained shallow (to B~ 25) photometry of individual stars in 19 clusters (of the 20 targets of the survey). We present the images 
and color magnitude diagrams (CMDs) of all of our targets. 

Methods. Point spread function fitting photometry of individual stars was obtained for all the WFPC2 images of the target clusters, 
and the completeness of the final samples was estimated using extensive sets of artificial stars experiments. The reddening, age, and 
metallicity of the clusters were estimated by comparing the observed CMDs and luminosity functions (LFs) with theoretical models. 
Stellar masses were estimated by comparison with theoretical models in the log(Age) vs. absolute integrated magnitude plane, using 
ages estimated from our CMDs and integrated J, H, K magnitudes from 2MASS-6X. 

Results. Nineteen of the twenty surveyed candidates were confirmed to be real star clusters, while one turned out to be a bright star. 
Three of the clusters were found not to be good YMC candidates from newly available integrated spectroscopy and were in fact found 
to be old from their CMD. Of the remaining sixteen clusters, fourteen have ages between 25 Myr and 280 Myr, two have older ages 
than 500 Myr (lower limits). By including ten other YMC with HST photometry from the literature, we assembled a sample of 25 
clusters younger than 1 Gyr, with mass ranging from 0.6 x 10 4 Mq to 6 x 10 4 Mq, with an average of ~ 3 x 1O 4 M0. Our estimates of 
ages and masses well agree with recent independent studies based on integrated spectra. 

Conclusions. The clusters considered here are confirmed to have masses significantly higher than Galactic open clusters (OC) in the 
same age range. Our analysis indicates that YMCs are relatively common in all the largest star-forming galaxies of the Local Group, 
while the lack of known YMC older than 20 Myr in the Milky Way may stem from selection effects. 

Key words. Galaxies: star clusters - Galaxies: individual: M31 - (Stars:) supergiants - Stars: evolution 
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1. Introduction 

Much of the star formation in the Milky Way is thought to have 
occurred within star clusters (Lada et al. I1991I Carpenter et al. 
12000b : therefore, understanding the formation and evolution of 
star clusters is an important piece of the galaxy formation puzzle. 
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Our understanding of the star cluster systems of spiral galax- 
ies largely comes from studies of the Milky Way. Star clusters 
in our Galaxy have traditionally been separated into two vari- 
eties, open and globular clusters (OCs and GCs hereafter). OCs 
are conventionally regarded as young (< 10 10 yr), low-mass 
(< 1O 4 M0), and metal-rich systems that reside in the Galactic 
disk. In contrast, GCs are characterized as old, massive systems. 
In the Milky Way, GCs can be broadly separated into two com- 
ponents: a metal-rich disk/bulge subpopulation, and a spatially 
extended, metal-poor halo subsystem (Kinman l 19591 Zinn[1985 
see also Brodie & S trader 2006] Harris|2001 for general reviews 
of GCs). 

However, the distinction between OCs and GCs has become 
increasingly blurred. For example, some OCs are luminous and 
old enought to be confused with GCs (e.g., Phelps & Schick 
2003). Similarly, some GCs are very low-luminosity systems 
(e.g., Koposov et al. 120071 1. and at least one has an age that is 
consistent with the OC age distribution (Palomar 1, Sarajedini 
et al. 120071 1. Moreover, a third category of star cluster, "young 
massive clusters" (YMCs) are observed to exist in both merg- 
ing (e.g., Whitmore & Schweizer 1995) and quiescent galaxies 
(Larsen & Richtler 1999). Indeed, YMCs have been known to 
exist in the Large Magellanic Cloud (LMC) for over half a cen- 
tury (Hodge fl 96 It . These objects are significantly more lumi- 
nous than OCs (My 5= -8 up to My ~ -15), making them 
promising candidate young GCs. Once thought to be absent in 
the Milky Way, recent observations suggest that their census may 
be quite incomplete, as some prominent cases have been found 
recently in the Galax y as well (Clark et al. 120051 Figer 120081 
Messineo et al. |20091 >. 

Thus, a picture has emerged that, rather than being distinct 
groups, OCs, YMCs and GCs may represent regions within a 
continuum of cluster properties dependent upon local galaxy 
conditions (Larsen 2003). The lifetime of a star cluster is depen- 
dent upon its mass and environment. Most low-mass star clusters 
in disks are rapidly disrupted via interactions with giant molec- 
ular clouds (Lamers & Gieles 2006 Gieles et al. 120071 1. These 
disrupted star clusters are thought to be the origin of much of 
the present field star populations (Lada & Lada 2003). Surviving 
disk clusters may then be regarded as OCs or YMCs, depend- 
ing upon their mass. Star clusters in the halo may survive longer 
since they are subjected to the more gradual dynamical processes 
of two-body relaxation and evaporation. The clusters which sur- 
vive for a Hubble time - more likely to occur away from the disk 
- are termed GCs (see also Krienke & Hodge 2007). To date, no 
known thin disk GCs have been identified in the Milky Way. 

After the Milky Way, M31 is the prime target for expand- 
ing our knowledge of cluster systems in spirals. However, our 
present state of knowledge about the M3 1 cluster system is far 
from complete. Similar to the Milky Way, M31 appears to have 
at least two GC subpopulations, a metal-rich, spatially concen- 
trated subpopulation of GCs and a more metal-poor, spatially 
extended GC subpopulation (Huchra et al. 1991 Barmby et al. 
2000). Also, again similar to the Milky Way GCs, the metal-rich 
GCs in M3 1 rotate and show "bulge-like" kinematics (Perrett et 
al. l2.002t . However, unlike the case in the Milky Way, the metal- 
poor GCs also show significant rotation (Huchra et al. 119911 
Perrett et al. |2002l Le eet al |2008l . Using the Perrett et al. ( [20021 
data, Morrison et al. (2004) identified what appeared to be a thin 
disk population of GCs, constituting some 27% of the Perrett et 
al. ((2002) sample. Subsequently, it has been shown that at least 
a subset of these objects are in fact young (< 1 Gyr), metal-rich 
star clusters rather than old "classical" GCs (Beasley et al. 2004 



Burstein et al. 2004 Fusi Pecci et al. [20031 Puzia et al.|2005 
Caldwell et al. 120091 . 

Fusi Pecci et al. (120051 hereafter F05) presented a compre- 
hensive study of bright young disk clusters in M31, selected 
from the Revised Bologna Catalog (RBC, Galleti et al. 120041 
by color [(B - V) < 0.45] or by the strength of the HB line in 
their spectra (HB > 3.5A). While these clusters have been noted 
since Vetesnik (1962) and have been studied by various authors, 
a systematic study was lacking. F05 found that these clusters, 
that they termed - to add to the growing menagerie of star cluster 
species - "blue luminous compact clusters" (BLCCs), are fairly 
numerous in M31 (15% of the whole GC sample), they have 
positions and kinematics typical of thin disk objects, and their 
colors and spectra strongly suggest that they have ages (signifi- 
cantly) less than 2 Gyr. 

Since they are quite bright (-6.5 ^ My ^ -10.0) and - at 
least in some cases - morphologically similar to old GCs (see 
Williams & HodgeEUOU hereafter WH01), BLCCs could be re- 
garded as YMCs, that is to say, candidate young GCs (see De 
Grijs 2009, for a recent review). In particular, F05 concluded 
that if most of the BLCCs have an age > 50 - 100 Myr they are 
likely brighter than Galactic open clusters (OC) of similar ages, 
thus they should belong to a class of objects that is not present, in 
large numbers, in our own Galaxy. Unfortunately, the accuracy 
in the age estimates obtained from the integrated properties of 
the clusters is not sufficient to determine their actual nature on 
an individual basis, i.e., to compare their total luminosity with 
the luminosity distribution of OCs of similar age (see Bellazzini 
et al. 2008 hereafter B08, and references therein). 

In addition to the question of the masses and ages of these 
BLCCs, it has become clear that the BLCC photometric and 
spectroscopic samples in M31 may suffer from significant con- 
tamination. Cohen, Matthews & Cameron (2006 hereafter C06) 
presented NIRC2@KeckII Laser Guide Star Adaptive Optics 
(LGSAO) images of six candidate BLCCs. Their K' very-high 
spatial resolution images revealed that in the fields of four can- 
didates there was no apparent cluster. This led C06 to the conclu- 
sion that some/many of the claimed BLCC may in fact be just as- 
terisms, i.e. chance groupings of stars in the dense disk of M3 1 . 
The use of the near infrared K' band (required by the LGSAO 
technique) may be largely insensitive to very young clusters that 
are dominated by relatively few hot stars, which emit most of 
their light in the blue region of the spectrum. Hence, the imaging 
by C06 may be inappropriate to detect such young clusters (see, 
for example, the detailed discussion by Caldwell et al. 2009). 
In any case, the study by C06 suggests that the true number of 
massive young clusters of M3 1 may have been overestimated. 

Therefore, in order to ascertain the real nature of these 
BLCCs we have performed a survey with the Hubble Space 
Telescope (HST) to image 20 BLCCs in the disk of M31 (pro- 
gram GO-10818, PL: J. Cohen). The key aims of the survey are: 

1 . to check if the imaged targets are real clusters or asterisms, 
and to determine the fraction of contamination of BLCCs by 
asterisms, 

2. to obtain an estimate of the age of each cluster in order 
to verify whether it is brighter than Galactic OCs of simi- 
lar age. Ultimately the survey aims to provide firm conclu- 
sions on the existence of a significant population of BLCCs 
(YMCs) in M31, in addition to OCs (see Krienke & Hodge 
120071 12008, and references therein) and GCs. 
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Fig. 1. Location of the 20 targets of our survey (empty circles) projected against the body of M31. The x symbols indicate the 
position of the additional ten Young Clusters we included in Sect. [4] 



In Perina et al. (2009a hereafter Pap-I) we have described in 
detail the observational material coming from our survey, and the 
data reduction and methods of analysis that we homogeneously 
adopt for the whole survey. We did that by taking the brightest of 
our surveyed clusters (VdBO) as an example. In this contribution 
we apply the same process to the whole sample, obtaining metal- 
licity, reddening and age estimates for all the targets of our sur- 
vey. We incremented our final sample of candidate M3 1 YMC by 
including in the final analysis ten further clusters having age es- 
timates available from the literature that are fully homogeneous 
with our own ones. In two companion papers, Hodge et al. (2009 
Pap-II, hereafter) identified and studied clusters of lower mass 



(with respect to those studied here) that were serendipitously im- 
aged in our survey, while Barmby et al. (2009, Pap-III, hereafter) 
studied the structure of the clusters that are the main targets of 
the survey. 

The paper is organized as follows. The sample is described 
in detail in Sect. 2, where we also summarize the data reduction 
procedure. In Sect. 3 we present the individual color magnitude 
diagrams (CMDs) and luminosity functions (LFs), we estimate 
ages, metallicities and reddening of each cluster. In Sect. 4 we 
derive the mass estimates for the clusters of our extended sam- 
ple (including data from the literature), we compare our clus- 
ters with open and globular clusters of the Milky Way and we 
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Fig. 2. F450W images of the 20 primary targets. Each image covers the central 10 x 10 on the PC field (10 = 38 pc at the 
assumed M31 distance modulus of 24.47). North is up and East to the left. 
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Table 1. Positional, photometric and spectroscopic parameters for the surveyed clusters. 



Name 


X" 
(arcmin) 


Y" 
(arcmin) 


R 

(arcmin) 


B 


V 




(B-V) (, ' WJ 


H/ U5 
(A) 


TT — /VH1 77? 

(A) 


B015D-D041 


-19.27 


9.22 


21.36 


19.11+0.02 


18.36±0.03 




0.15 


7.32 


l 


B040-G102 


-35.40 


-11.92 


37.35 


17.54±0.03 


17.20±0.04 


0.18 


0.11 


7.41 


7.58± 0.30 1 


B043-G106 


-33.62 


-11.37 


35.49 


17.04±0.03 


16.77±0.04 


0.17 


0.04 


5.53 


5.70± 0.30 1 


B066-G128 


-29.55 


-13.17 


32.35 


17.56±0.03 


17.35±0.04 


0.25 


-0.02 


4.67 


4.84± 0.30 1 


B081-G142 


-25.26 


-12.36 


28.12 


17.36±0.02 


16.86±0.03 


0.43 


0.20 


7.98 


8.15±0.30 1 


B257D-D073 


45.98 


4.02 


46.16 


18.41±0.02 


18.00±0.04 




0.01 


5.49 


5.66± 0.30 1 


B318-G042 


-52.14 


-1.32 


52.16 


17.02±0.03 


16.82±0.03 


0.06 


0.03 




5.49±0.12 1 


B321-G046 


-55.50 


-7.41 


55.99 


17.82±0.02 


17.51±0.03 


0.11 


0.06 


6.29 


6.85+ 0.32 1 


B327-G053 


-47.67 


-3.45 


47.79 


16.75±0.03 


16.58±0.03 


0.21 


-0.03 


4.09 


3.78±0.14 1 


B376-G309 


42.16 


-10.67 


43.49 


18.35±0.02 


17.97±0.04 


0.34 


0.08 




6.40± 0.06 1 


B448-D035 


-43.16 


-2.97 


43.26 


18.01±0.03 


17.46±0.04 


0.50 


0.20 


6.70 


6.87± 0.30 1 


B475-V128 


45.00 


4.06 


45.18 


17.55±0.03 


17.09±0.04 


0.20 


0.11 


5.96 


6.13±0.30 1 


Viol 




7 1 1 
1 . 1 / 




i q i ^_i_n ni 
lo. 1D±U.Uj 


1 1 AOj.fl (\A 


U.j / 




J. (54 


c. Al i r\ on i 
O.U1+ U.jU 1 


B083-G146 


19.83 


22.08 


29.68 




17.09'' 


0.65 


0.56 


3.75 


1.75± 0.42 1 


B222-G277 


10.22 


-16.16 


19.12 


18.00±0.02 


17.24±0.03 


0.57 


0.56 


8.47 


4.46+0.31 1 


B347-G154 


27.74 


26.74 


38.53 


17.23 d 


16.50'' 


0.62 


0.67 




2.87+0.17 2 


B374-G306 


41.13 


-10.55 


42.46 


18.69±0.03 


18.23±0.04 


0.33 


0.16 


4.07 


4.24+ 0.30 1 


NB16 


1.96 


4.19 


4.63 


18.83+0.04 


17.59±0.10 


0.55 


0.99 




3.34+ 0.08 2 


VDBO 


-47.16 


-4.33 


47.36 


14.94+0.09'' 


14.67±0.05 c 


0.12 


0.07 


4.30 


4.50+ 0.07 1 


NB67-AU13 


1.68 


3.73 


4.09 


16.48±0.02 


15.92±0.03 


0.37 


0.36 







B and V magnitudes are from new aperture photometry performed on the CCD images of Massey et al. (2006), except for B083 and B347 that 

are not included in the area covered by that survey. 
a X and Y are projected coordinates in the direction along (increasing Eastward) and perpendicular to the major axis of M31, in arcmin. 
b ff is a flag indicating if the target has been selected from Table 1 or Table 2 of F05. 
c From Pap-I. 
d From the RBC. 

w '' from this work: B and V from this table and E(B-V) as estimated in Sect. 3 from isochrone fitting. 

15 from Fusi Pecci et al. (2005): (B-V) are calculated assuming a single value of E(B-V)=0.1 1 for all the clusters. 

59 from Galled et al. (2009). 



compare our estimates with those from the recent and exten- 
sive analysis of young M31 clusters by Caldwell et al. (2009 
hereafter C09), that are based on integrated spectra. In Sect. 5 
our main results are briefly summarized and discussed. Finally, 
in Appendix A we report on M3 1 clusters or candidate clusters 
listed in the RBC that have been serendipitously imaged within 
our survey, and, in Appendix B, we report on the nature of candi- 
date BLCC=YMC M31 clusters that have an HST image in the 
archive, independent of this survey. 

2. Description of the sample 

Table [TJ lists the target clusters of our survey and reports some 
positional and spectro-photometric parameters that were rele- 
vant for their selection. New homogeneous large-aperture (r ap ~ 
5" - 10", depending on the curve of growth of each cluster) 
integrated B,V photometry for all the targets has been obtained 
from the publicly available CCD images by Massey et al. (2006), 
and calibrated using the published photometry from the same au- 
thors, as done in Pap-I for VdB-0 (see Pap-I for further details). 

Fig-fflshows that the vast majority of the targets are projected 
onto the so-called 10 kpc ring (see Hodge 1992 Barmby et al. 
2006 , C09 and references therein), a site of ongoing star forma- 
tion in the thin disk of M3 1 . The only exceptions are B347 and 
B083, that are significantly farther from the center of the galaxy, 
and NB16 that is projected onto the outer regions of the M31 
bulge. We will see below that these three clusters do not fulfill 
the selection criteria by F05 for bona fide candidate YMCs and, 
in fact, they are likely old (see Sect. 13.3b . 



Eighteen of the twenty targets were drawn from Tab. 1 of 
F05, i.e. they were confirmed cluster^] that were classified as 
genuine BLCC = YMC by these authors as they had Hp > 3. 5 A 
or, when lacking a measure of Hp, (B - V)o < 0.45. After a care- 
ful inspection of the HST archive, we excluded from the selec- 
tion any cluster from Tab. 1 of F05 that had already been imaged 
with HST (serendipitously, in most cases, see Appendix B), and 
we chose the brightest 18 among the remaining ones. F05 as- 
sumed E(B-V)= 0.11 for all the considered sample, in Sect. 3 
we will show that the typical reddening of these clusters is sig- 
nificantly higher than this, in most cases E(B-V)> 0.20, in good 
agreement with the estimates by C09 (see Fig. [17}. Hence, in 
general, the (B - V)o colors derived here are bluer than those 
adopted by F05. Galleti et al. (2009, G09 hereafter) presented 
new estimates of the Hp index (with respect to those reported 
by F05), taken either from their own observations or from the 
recent literature. In Table [TJ we report both the (B - V)o and Hp 
values from F05 (that were used for the selection of the sam- 
ple) and those derived here and in G09, when availablfl In 
one case (B083) the new value of Hp is much lower than that 
reported by F05 (1.75A instead of 3.75A) and than the selec- 
tion limit. Moreover, even with the new E(B-V) estimate derived 



2 RBC class f=l, meaning that they have been classified as bona- 
fide M31 clusters by some author, based on their spectra and/or high 
resolution images. 

3 Note that the scales of the Hp index adopted by F05 and G09 are 
slightly different. The Hp > 3.5A threshold by F05 translated into Hp > 
3.1 A in the scale by G09 (see the latter paper for discussion and details). 
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here, (B — V)o — 0.551, significantly redder that the limit adopted 
for the selection. For these reasons B083 can no longer be con- 
sidered as a candidate YMC, as it does not fulfill the selection 
criteria when the newly available data are considered. The analy- 
sis of the CMD (in Sect. 3) will confirm that the cluster is in fact 
much older than genuine YMC, and possibly as old as classical 
GCs. 

The remaining two targets (NB16 and B347) were selected 
form Tab. 2 of F05, including clusters not fulfilling their selec- 
tion criteria for YMC but classified as young (or possibly young) 
by some author in the past. In both cases Hp were lacking at the 
time, and the new values reported by G09 are significantly be- 
low the selection threshold for a YMC. B347 is also much redder 
than (B-V)o = 0.45. On the other hand, we find (B-V) = 0.399 
for NB16. In this case the criterion based on Hp must prevail 
over that based on de-reddened color as the former is reddening- 
independent, while relatively low photometric and/or reddening 
errors can shift the color of this cluster above or below the se- 
lection threshold. In conclusion, the newly available data indi- 
cates that both NB16 and B347 are not good YMC candidates, 
as will be confirmed by their CMDs (see Fig.[T2b. Hence, just re- 
considering the original selection in the light of new estimates of 
integrated properties, our sample of bona fide YMC candidates 
is reduced to 17 objects, including VdBO which was studied in 
detail in Pap I. 

Postage stamp images of all the targets, from our HST data, 
are presented in Fig. [2] (see Sect. 12. 11 1. Inspection of the images 
reveal that all our targets are actually genuine clusters, with the 
only exception of NB67 that is a bright star projected into a 
dense background of M3 1 (disc) stars (see also Pap-III, for the 
light profiles of the clusters). For obvious reasons NB67 will 
be not considered further in the following analysis. A first con- 
clusion that can be drawn just from this preliminary analysis is 
that the incidence of spurious objects in our sample is of 1/17^ 
6%, much lower than hypothesized by C06. If we consider the 
set of 36 objects listed by F05 in their Tab. 1 for which HST 
images were available in the archive we obtain the same result 
(see Appendix iBl for discussion and further details). Moreover, 
none of the considered clusters is in fact an asterism (includ- 
ing those considered in Appendix iBlPl Finally, if we extend our 
analysis to all the objects classified as YMC by F05 that have 
been ever imaged with HST we find the same very low degree 
of contamination (see Appendix IBl. Hence we are dealing with 
a significant class of real stellar systems. A second conclusion is 
that while some of the considered cluster appear quite extended 
and sparse (like, for example, B257D, B475, and V031), there 
are also rather compact globular-like clusters (like, B043, B081, 
and B327, as noted earlier B347 is likely old). 

2.1. Observations, data reduction and assumptions. 

The characteristics of the survey data and the whole process of 
data reduction and data analysis that has been applied in this 
study is described in detail in Pap-I. In these section we briefly 
summarize the key characteristics of the dataset and of the pro- 
cess, for the convenience of the reader. 

Two t exp = 400 s images per filter (F450W and F814W) 
were acquired for each cluster with the Wide Field and Planetary 
Camera (WFPC2) on board of HST, keeping the target at the 
center of the PC field. Unlike the case of VdBO, treated in Pap-I, 
the clusters studied here have limiting radii significantly smaller 



4 Bright stars are well-known classical contaminants in lists of can- 
didate M3 1 clusters of any kind, see Galled etal. l2006al 




F814W 



Fig. 3. Completeness (Cf) of the samples as a function of 
F814W magnitude, obtained from artificial stars experiments, 
for all the clusters of our survey (listed in Tab. table: 1) and for 
two different color ranges. The upper panel is for a color range 
enclosing the MS of young clusters, the lower panel is for a color 
range enclosing the red giant stars. The C /(F814W) function of 
each cluster (for each color range) is computed considering only 
artificial stars enclosed in the radial range that is used to select 
the sample dominated by cluster stars that will be studied in the 
following (typically r < 5", see Sect. l2.2l and Sect. 01. Note that 
all the C f(F8 14W) functions are very similar, except for the case 
of the exceedingly compact (and crowded) cluster NB 16, labeled 
in both panels. 



than the size of the PC camera (^ 39" x 39", see Pap-III), there- 
fore both the cluster population and the surrounding field can be 
studied using the PC images alone (see Sect. 12.2b without rely- 
ing on the WF cameras. The analysis of the field population in 
the portions of the M3 1 disk sampled by our WF images will be 
the subject of another contribution (Perina et al., in preparation). 

Photometry of the individual stars has been obtained with 
HSTPHOT (Dolphin I2000ab . a Point Spread Function fitting 
package specifically developed for WFPC2 data. The reduc- 
tion process includes cleaning of cosmic -ray hits and bad pix- 
els, correction for Charge Transfer Efficiency (CTE, Dolphin 
2000b), and absolute photometric calibration in the VEGAMAG 
system (Holtzman et al. 119951 Dolphin I2000bl >. The images 
were searched for sources having peak intensities at 3<x above 
the background. The output catalogs were cleaned of spuri- 
ous and/or badly measured sources by selecting stars with 
HSTPHOT global quality flag=l, crowding parameter < 0.3, 
X 2 < 2.0 and \sharp\ < 0.5. The final catalogs containing po- 
sition and F450W, F814W photometry of the PC fields will be 
made publicly available through a dedicated WEB page^. 

We estimated the completeness of our samples as a function 
of magnitude, color and position on the field by means of ex- 
tensive artificial stars experiments (more than 10 5 artificial stars 



5 www.bo.astro.it/M31/YMC 
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Fig. 4. Selection boxes used for the stellar surface density pro- 
files shown in Fig. [5] and [6l are superimposed on the CMD of 
two of the surveyed clusters taken as examples: a young clus- 
ter with a prominent MS (left panel) and an older cluster dis- 
playing just the tip of the RGB (right panel). The blue box at 
F450W - F814W ~ 0.5 selects bright MS stars (young popula- 
tion), the faint redder box (F450W - F814W > 1.0) selects red 
giant stars (old population). In a few cases, the boxes have been 
slightly shifted in color to best match the MS and RGB features 
of a cluster with higher reddening. 



were simulated, per field of view, i.e. more than 4 x 10 5 per clus- 
ter), as described in detail in Pap-I. Fig.[3]show the completeness 
factor (Cf) as a function of magnitude for all the clusters, for 
two different color ranges (one covering the clusters' main se- 
quence (MS) and one covering the Red (Super) Giant branches). 
The reported Cf curves refers to the circles enclosing most of 
the cluster population that are defined in Sect. 12.21 hence they 
are fully relevant for the following analysis. Note that the com- 
pleteness conditions are very similar for all the clusters (includ- 
ing VdBO, presented in Pap-I), except NB16. This cluster is so 
compact that the considered region is much more crowded than 
all the other cases, thus the completeness is significantly worse. 
The typical photometric uncertainties as derived from the artifi- 
cial stars experiments are < +0.02 for F450W F814W < 21, 
< ±0.05 for F45QW F814W < 22.5, and ±0.2 for 
F450W - F814W < 24.0 (see Pap-I, for details). 

In the following we will always assume (m - M)o = 24.47, 
from McConnachie et al. (2005), corresponding to D = 783 kpc. 
At this distance 1" corresponds to 3.8 pc, 1' to 228 pc. We adopt 
Aptfow = 4.01 5E(ff - V) and A FUm = 1.948£(fi - V), from 
Schlegel et al. J 1998b . We will use theoretical isochrones and 
LFs in the HST/WFPC2 VEGAMAG system from the set by 
Girardi et al. (2002, hereafter G02), considering only models in 
the range of metallicity |Zq ^ Z ^ 2Zq, that seem appropri- 
ate for young disk clusters. Details and discussion regarding the 
choices outlined above can be found in Pap-I. 



2.2. Radial selection and first classification 

Before proceeding with the analysis of the CMDs of the clus- 
ters, we need to select - for each cluster - a sub-sample of the PC 
field that is as representative as possible of the cluster popula- 
tion, possibly minimizing the contamination by the surrounding 
M31 field. Following Pap-I we adopt a radial selection, retaining 
in the final cluster sample the stars lying within a certain distance 
from the cluster center. To determine the selection radius to be 
adopted for each individual cluster we proceeded as follows: 

- We defined two broad selection boxes on the CMD, one en- 
closing the bright MS typical of young clusters (Blue Box) 
and one enclosing a redder region that should be dominated 
by old stars at the tip of the red giant branch (RGB) but 
can enclose also intermediate-age asymptotic giant branch 
(AGB) and some red super giant (RSG) stars, as illustrated 
in Fig. |4] (Red Box). 

- We derived surface-density radial profiles by counting stars 
selected in the two boxes on concentric annuli. To obtain 
smoother profiles with the relatively low number of stars 
available we adopted overlapping annuli of width 1.8", with 
a radial step of 0.9" between subsequent annuli. The profiles 
from main sequence (MS) stars and from red stars (shown in 
Fig|5]and|6]l are normalized to the minimum surface-density 
encountered in the raster of radial annuli, that should be con- 
sidered as roughly representative of the surrounding field. 
For example, the profiles of B066, in the middle left panel 
of Fig [5] shows that at the center of this cluster the surface 
density of bright MS stars is ^ 20 times higher than in the 
surrounding field, while there is no overdensity of red stars 
correlated to the cluster. 

- Based on the scale of the detected overdensity we fixed the 
selection radius of each cluster (marked in the plots as a 
vertical dashed line), with the aim of isolating a circle that 
should be dominated by cluster stars. The typical selection 
radius is r ~ 5". 

In the following we will analyze only the CMDs of the radi- 
ally selected samples, as the best representation of the population 
of each cluster. The CMDs of the surrounding fields are shown 
in Fig. [7] for comparison with those of the respective clusters 
that are studied in detail in Sect. [3] 

Fig [5] and [6] deserve some further comment. First of all, it 
has to be noted that all the clusters (at their centers) show an 
overdensity of a factor of £ 10 with respect to the surrounding 
field, at least in one of the two profiles. The only exception is 
NB16 that is so compact that only a tiny corona is resolved into 
stars, resulting in a low (~ 2x) overdensity of red stars (but see 
the light profile obtained in Pap-III). Note that in many cases, the 
very central region of the cluster is not fully resolved, thus the re- 
ported central overdensities are just lower limits to the true ones. 
Second, there are five clusters that show no sign of overdensity 
in the Blue Box. B083, B347, and NB16 have been discussed 
above; they cannot be considered as YMC candidates anymore. 
B222 and B374 on the other hand have both H p > 3. 5 A. In four 
cases the cluster show no sign of overdensity in the Red Box, in 
particular, B040, B043, B066, B327. In all the other cases, the 
overdensity is detected in both the Blue and Red boxes popula- 
tions, even if not necessarily in similar degree. In general the 
overdensity from MS stars is larger than in RGB/AGB/RSG, 
as expected from evolutionary considerations (Renzini & Fusi 
Pecci 1981. 
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Fig. 5. Stellar surface density profiles of the young (open circles connected by a continuous line) and old (crosses connected by a 
dashed line) populations (as defined by the selection boxes illustrated in Fig. |4j for nine of the surveyed clusters. 



3. Age and metallicity 

Once established that our targets are real clusters, the main pur- 
pose of our survey is to obtain a reliable age estimate for all 
of them from their CMDs. This will be done by comparison 
with theoretical isochrones from the set by Girardi et al. (2002 
G02 hereafter, the models are in the same photometric system 
as the data; see Pap-I for a discussion about the choice of the 
set of theoretical models), following the approach described in 
detail in Pap-I. The procedure provides a simultaneous estimate 
of the age, the reddening and the metallicity of each cluster un- 
der consideration, by eye-aided isochrone fitting. In Pap-I we 



have shown that the data from our survey can be used to reli- 
ably estimate ages in the range from ~ 10 Myr to < 500 Myr 
(also depending on the total mass of the considered clusters, i.e. 
on the number of stars populating the MS), from the luminosity 
and color of the Turn Off (TO) point. The distribution of RSG 
may help to constrain the metallicity of the population, while 
the color of the blue edge of the MS is the best indicator of the 
degree of interstellar extinction (see Pap-I). 

In our sample, there are eleven clusters that have a signifi- 
cant number of MS stars brighter than FS14W = 24.0. As the 
completeness of the sample is Cf ~Z 80% above this limit, (in the 
color range enclosing the MS, see Fig. [3]), reliable completeness- 
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Fig. 6. Same as Fig.[5]for the remaining nine surveyed clusters. 



corrected LFs of the MS population can be obtained, and used 
to further constrain the age of these clusters, as one in Pap-I. All 
of these eleven clusters have ages lower than 200 Myr. They 
are homogeneously analyzed in Sect. 13.11 Also VdBO belongs 
to this class but it is not considered here as it has been already 
treated in Pap-I. 



Finally, there are five clusters that do not display any obvious 
MS population in the range of magnitudes accessible with our 
data. For these clusters we can provide only a strong lower limit 
to their age, that must be older than 300-500 Myr. These clusters 
are discussed in Sect. 13.31 The final results of the analysis of the 
CMD presented below are reported in Tab. [2] 



Two clusters (B475 and V031) show a clear MS population 
only for F814W > 24.0. As their observed MS lie in a range 
where the completeness factor drops from Cj ~ 80% to Cf ~ 
in ~ 2 magnitudes their LF would be strongly affected by large 
completeness corrections. For these reason we limit our analysis 
to isochrone fitting for these clusters (Sect. [3~!2"l ). 



3.1. Clusters with bright MS (age< 200 Myr) 

Fig. [8] |9] and [10] show the observed CMDs and LFs of the 
eleven clusters having a significant MS population brighter than 
FS14W = 24.0. The boxes overplotted on the CMDs have been 
used to select the stars that were used to derive the LFs. 
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Fig. 7. CMDs of the fields surrounding the target clusters. Only stars lying in the radial range 5"r < 16.5" on the PC chips are 
plotted. The thin lines are the loci where the completeness reaches 50%. 



For each cluster we explored the space of parameters to find 
the isochrone and the reddening providing the best overall fit to 
the observed CMDs. As differential reddening may move stars 
toward the red and the presence of binary systems also has the 
effect of broadening the MS toward the red side, we searched for 
solutions where the theoretical MS fits the blue side of the MS. 



As noted above, the distribution of RSGs was used as a guide 
to fix the metallicity of the best-fit model (see Pap-I). Following 
the approach of Pap-I, we adopt Z=0.019 as the starting guess for 
the metallicity of the cluster, trying other metallicity only if this 
was required to better fit some feature of the CMD. A correct in- 
terpretation of the cluster CMD was aided by a comparison with 
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Fig. 8. Left panels: CMDs of the clusters B327, B015D, B066, and B318, displaying only stars within the radial selection reported in 
the upper right corner of each panel. The adopted best-fit value of the reddening and the age and metallicity of the best-fit isochrone 
(thick continuous line) are reported in the lower right corner of each panel. The rectangular boxes adopted to select the stars used to 
obtain the LFs shown in the right panels are also plotted. Right panels: the observed completeness-corrected LFs of the cluster MS 
(filled circles with error bars) are compared with theoretical models of different ages. The thick continuous line corresponds to the 
best-fit model shown in the CDMs. In all cases, it provides a reasonable fit to the observed LF and, in particular, to the sudden drop 
of star counts at the upper limit of the MS. The dotted and dashed lines are theoretical LFs corresponding to strong upper and lower 
limits to the age, respectively, as they are the nearest models that can be clearly excluded by the data. The theoretical LFs have been 
arbitrarily normalized to best match the three faintest observed points. 



the CMD of the surrounding field, to establish, for example, if a 
population of a few RSG can be considered as characteristic of 
the cluster or compatible with belonging to the field. The typical 
uncertainty on the reddening estimate is +0.04 mag (see Pap-I). 



The theoretical LF of the isochrone that best-fits the ob- 
served CMD morphology (thick continuous line in the right pan- 
els) is compared to the observed LF (filled dots with error bars) 
to check the compatibility of the solution with the star counts 
(Salpeter's 119551 Initial Mass Function is adopted). In all the 
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Fig. 9. Same as Fig.|8]but for the clusters B040, B043, B257D, and B448. 



cases considered the adopted theoretical LF is in good agree- 
ment with the observations and, in particular, it reproduces the 
sudden drop in star counts corresponding to the upper luminos- 
ity limit of the MS, a feature that is mainly sensitive to age (see 
Pap-I and references therein). Two theoretical LFs of the same 
metallicity as the main solution but different ages are used to 
show the maximum and minimum age that are not compatible 
with the observed LF. The difference between these values and 
the age of the best-fit solution are taken as the uncertainty as- 
sociated with our age estimate. Nine of the eleven clusters con- 
sidered in this section have ages between 50 Myr and 100 Myr. 
All of them show a recognizable (and in same case sizable, see 
B040, for example) population of RSG stars, in addition to an 



obvious MS. The other two clusters, B081 and B321 have ages 
of 140 and 170 Myr, respectively. 

3.2. Clusters with faint MS (200 Myr<age< 500 Myr) 

Fig. QT| shows the CMDs of the two clusters whose MS is 
fainter than F814W = 24.0. The F450W magnitude is plotted 
here instead of F814W (adopted in Fig. |8] |2| and [TOj as this 
makes the faint MS of these clusters more clearly visible. The 
best fit isochrones are plotted as thick lines. The thin lines are 
isochrones having ages that bracket the age solutions that can be 
considered still compatible with the data. The difference in age 
between these solutions and the assumed best-fit are adopted as 
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the uncertainty associated with our age estimates for this cases 
(see Pap-I). The two clusters have ages of ^200 Myr (B475) and 
=280Myr(V031). 



3.3. Clusters whose MS is not detected (age> 500 Myr) 

Fig. [12] shows the CMDs of the clusters that do not display a 
clear MS in the considered range of magnitudes. In each panel 
we plot (a) the "youngest" isochrone that is compatible with the 
observed CMD morphology, to provide a firm lower limit to the 
age of these clusters (thick continuous line), and, (b) a 12 Gyr old 
isochrone (thick dashed line), showing that the observed CMD is 
also compatible with very old ages. In all the cases we adopt the 
metallicity value that provided a satisfactory match of the color 
of the (putative) RGB. 

Three of the five clusters considered here (B083, NB16 and 
B347) have integrated properties that are compatible with old 
ages (see Sect.|2]l. B083 and B347 display a steep and well pop- 
ulated red sequence, much bluer than the limits imposed by the 
run of the completeness as a function of color (thin dotted lines), 
typical of the RGB of classical old (and metal deficient) GCs. 
The handful of stars resolved in NB16 are also compatible with 
being near the tip of an old RGB, but their scarcity poses strong 
caveats on any interpretation. 



B347 and B222 are more interesting cases: both have two 
independent concordant estimates of Hp indicating Hp > 4.0A, 
and both have some stars just above the detection limits in the 
blue, that may be compatible with the bright end of a fainter 
MS. The observational scenario is fully consistent with the hy- 
pothesis that these two clusters might be intermediate-age (age~ 
0.5 - 2 Gyr). A deeper photometry follow-up is clearly required 
to settle the issue of the age of these clusters. It is worth noting 
that a convincing case for an M31 cluster in the age range 1-8 
Gyr with age estimated from a CMD has never been provided. 

4. Masses from ages and J,H,K integrated 
photometry 

In Table|2]we report the age, metallicity and reddening estimates 
obtained from the analysis of the CMDs presented above. To in- 
crease the sample of YMC to be considered in the following 
we added a total of 10 further clusters whose ages have been 
derived from CMDs obtained from HST data in a way fully ho- 
mogeneous with that adopted here. In particular we add six clus- 
ters from Perina et al. (2009b, P09b hereafter) and four clus- 
ters from Williams & Hodge (2001, WH01 hereafter; see Pap-I). 
All of them lie in the range of V luminosities typical of YMC 
(My ^ -6.5, according to F05), with the only (possible) excep- 
tions of M050 and M039 that appear somewhat fainter than this, 
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Fig. 11. Observed CMDs of the clusters B475 (left panel) and V031 (right panel) in the plane F450W vs. F450W-F814W where the 
MS population of these older clusters is more clearly visible. Only stars with the radial selection reported in each panel are plotted. 
The best-fit isochrone is plotted as thick line (age, metallicity and reddening values are reported in each panel). The thin isochrones 
bracket the upper and lower limits on the age, and correspond to age 125 Myr and 315 Myr for B475, and age 200 Myr and 400 
MyrforV031. 



and of B521 that lacks an estimate of its V magnitude (but it is 
found to have a mass similar to other YMC, based on its Near 
Infrared Magnitudes, see below). We decided to keep these clus- 
ters within our sample, being well aware that the threshold be- 
tween the brightest of the clusters studied in Pap-II and Krienke 
& Hodge ( 2007 , 2008) and the faintest clusters considered here 
is somewhat blurred, both by lack of a clear-cut definition and 
by observational uncertainties. In particular, Fig. [20] will show 
that some of the clusters studied in Pap-II appear to have masses 
typical of YMC. Still we preferred not to include these massive 
Pap-II clusters as main objects of the present analysis as most of 
them have their ages estimated from integrated colors, i.e. with 
significantly greater uncertainties than those obtained here from 
CMDs (see, e.g., Fig. 8 of Pap-Ilfl 



6 There are only two clusters from Pap-II having M v S -6.5 and 
ages estimated from their CMD, but also in these cases the associated 



Five of the newly included clusters are projected onto the 10 
kpc ring, as most of our original targets, four lie slightly nearer to 
the center of the galaxy, and one is in the outskirts of the visible 
disk (see Fig.Q). B049, B367, B458, B315 and B317 have two 
independent estimates of Hp, all of them higher than 4.5 A (F05, 
G09). B342 has just one estimate (Hp = 7.06A, FP05), while the 
other four clusters lack any measure of this index. B368 lacks 
Hp but has (B - V) = 0.06. For M039, M050 and B521 there 
is no (B - V)o estimate available. In any case all the six clusters 
from P09b and the four from WH01 have age < 1 Gyr, as derived 
from their CMD. 

To derive the most reliable estimate of the total stellar mass 
of the clusters in our sample we couple our age estimates with 
integrated Near Infra Red (NIR) photometry, as stellar mass-to- 
light ratios in NIR bands have a much shallower dependence on 



age uncertainties are relatively large, i.e. 0.5-0.6 dex in log(Age) vs. a 
typical uncertainty of 0.2 dex for our main sample, see Tab. [2] 



S. Perina et al.: A HST/WFPC2 survey of bright young clusters in M31. II. 



15 



'I 1 1 1 I 1 

B374 

E(B-V)=0.30 
Z=0.019 
t=300 Myr 



r<3" 



1 I 1 1 1 I 
B222 

E(B-V) = 0.20 
Z=0.019 
t=400 Myr 



I 1 1 1 I 
B083 

E(B-V)=0.20 
Z = 0.008 
t=500 Myr 



20 - 



2 22 
oo 



24 




I | I ' I I | I I I 



r<4" 





20 



- 22 
oo 



NB16 

E(B-V)=0.25 
Z=0.019 
t=500 Myr 



r<1.5" 



B347 

E(B-V) = 0.06 
Z=0.008 
t=630 Myr 



r<6" 



2 A 

F450W-F814W 



24 - 




2 4 

F450W-F814W 



2 4 

F450W-F814W 



Fig. 12. CMDs of the clusters B374, B222, B083, NB16, and B347. Only stars within the radial selection reported in each panel 
are plotted. The thin dashed lines marks the locus where the completeness of the sample reaches 0% (see Pap 1), to illustrate 
the selection effects on the CMD morphology imposed by the run of limiting magnitude as a function of color. In each panel, the 
continuous line is the youngest age isochrone that is compatible with the observed CMD, providing a strong lower limit to the age 
of each cluster. The adopted age, metallicity and reddening values are reported in the upper left corner. The dashed line is a 12 
Gyr old isochrone matching the color of the observed RGB. The metallicity of these old-age isochrones is Z = 0.001, 0.004, 0.001, 
0.004, and 0.001 for B374, B222, B083, NB16, and B347, respectively. 



age than their optical counterparts (see Pap-I for discussion). As 
the best estimate of the integrated J,H,K magnitudes we took the 
values of the r = 10" aperture magnitudes from the 2MASS- 
6X-PSC catalog (see Nantais et al. |20061 l, that is obtained from 
deeper observations (with respect to the normal 2MASS data, 
Skrutskie et al. 120061 1 over a limited region of the sky that, 
luckily, includes M3 1 . The adopted NIR photometry as well as 
the accurate positions reported in 2MASS-6X-PSC are listed in 
Table [3] Only two clusters have no valid measures in 2MASS- 
6X-PSC, i.e. B367 and M039. To preserve the homogeneity of 
the analysis we do not include these clusters in any of the follow- 
ing analyses that make use of mass estimates, however, for com- 
pleteness, in Tab.[3]we provide a tentative mass estimate derived 
from the log(age) vs. My diagram presented in Fig. [14] The ap- 
parent magnitudes are transformed into absolute ones adopting 
the reddening estimates derived here (Tab. the distance mod- 



ulus (from McConnachie et al. 120051 and the reddening laws 
(from Rieke & Lebofsky 1985) adopted in Pap-I. 

In Fig. [13] we compare the position of our clusters in the 
integrated (J,H,K) magnitude vs. log(age) plane with a grid of 
models of Simple Stellar Population (SSP) of solar metallicity 
and various total mass, from the set by Maraston (1998 2005 
see Pap-I). In B08 and in Pap-I we have shown that the mass that 
can be deduced from these plots depends only weakly on the as- 
sumed metallicity and IMF. Here we get an independent estimate 
of the mass from each (J,H,K) plot and we take the weighted av- 
erage of the three values as our final estimate. The uncertainties 
were obtained on each individual estimate from J, H, K by find- 
ing the maximum interval in mass that was compatible with the 
errors in age and in integrated magnitudes. Then the three values 
(per cluster) were combined into the final weighted error that is 
reported in Table [3]together with the final mass estimates. 
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Table 2. Newly derived ages, metallicity and reddening for the target clusters and other clusters included in the analysis". 



Name 


log(t) 


Alog(t) 


z 


E(B-V) 


M v * 


This survey 


B015D-D041 


7.85 


±0.15 


0.019 


0.60 


-8.53 


B040-G102 


7.90 


+0.20 
-0.15 


0.019 


0.23 


-7.80 


B043-G106 


7.90 


+0.20 


0.019 


0.23 


-8.22 


B066-G128 


7.85 


±0.15 


0.019 


0.23 


-7.76 


B081-G142 


8.15 


±0.15 


0.019 


0.30 


-8.60 


B257D-D073 


7.90 


+0.20 


0.019 


0.40 


-8.31 


B318-G042 


7.85 


±0.15 


0.008 


0.17 


-7.98 


B321-G046 


8.23 


+0.10 
—0.15 


0.019 


0.25 


-7.57 


B327-G053 


7.70 


+0.15 
-0.10 


0.008 


0.20 


-8.51 


B376-G309 


8.00 


±0.15 


0.019 


0.30 


-7.34 


B448-D035 


7.90 


+0.20 
-0.15 


0.019 


0.35 


-8.07 


B475-V128 


8.30 


±0.20 


0.008 


0.35 


-8.00 


V031 


8.45 


±0.15 


0.004 


0.35 


-8.12 


VDBO 


7.40 


±0.30 


0.019 


0.20 


-10.03 


B083-G146 
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For five surveyed clusters only a lower limit to the age can be obtained from our CMDs. 
a The additional clusters are six clusters studied in Perina et al (2009a), from HST archive data, and the four clusters studied by Williams & 

Hodge (2001). 
b Integrated V magnitudes from the RBC. 



It is very reassuring to note that the three plots provide very 
similar age estimates: all the clusters considered appear to have 
masses between ~ 1O 4 M0 and ~ 1O 5 M0. The estimates from 
the three different NIR magnitudes typically agree within a fac- 
tor of 2. The adoption of a Kroupa (2001 ) IMF instead of that of 
Salpeter would change the mass estimates by less than a factor 
of 2 (Pap-I). The adoption of different sets of models would lead 
to a maximum difference of the same amount in the final mass 
estimates (we have compared the M/L predictions adopted here 
with those from the sets by Pietrinferni et al. 120041 and Bruzual 



& Chariot 2003 , in the age range that is relevant for our clus- 
ters). Finally, if models with age-dependent M/L are adopted 
(i.e. including the effects of differential mass loss, Kruijissen& 
Lamers 2008 ), the mass estimates for our clusters change by a 
mere ^ 20% (see also Pap-III). Taking all of these factors into 
account it turns out that our mass estimates should be accurate 
within a factor of 5 3, as confirmed also by the comparison with 
the independent estimates from Pap-III and C09. 

There is only one case of significant disagreement in the 
position of a cluster in the different NIR passbands, i.e. B347 
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Fig. 13. Log(age) vs. integrated magnitude plane for near infrared colors. The target clusters are represented as open squares (VDBO 
as a crossed square), the clusters from P09b as open stars, and the clusters from WH01 clusters as open triangles, IR magnitudes 
are taken from Tab. [3] Note that B367 and M039 are not plotted because they lack NIR photometry in the 2MASS-6X-PSC catalog. 
The gray symbols show the clusters that have "null" error on IR magnitudes in the 2MASS-6X-PSC catalog. Integrated magnitudes 
of Galactic GCs (x symbols) are taken from Cohen et al. (2007). The continuous lines are fixed-stellar-mass models from the set 
by Maraston (1998, 2005) for SSPs of solar metallicity, with a Salpeter's Initial Mass Function (IMF) and intermediate Horizontal 
Branch morphology. Note that in this plane, the dependence of the models from the assumed IMF, metallicity and HB morphology 
is quite small (see B08). The dotted lines are M — 10 4 Mq and M - 10 5 Mq iso-mass models assuming a Kroupa 2001 IMF instead 
of a Salpeter ( 1955) IMF, plotted here to illustrate the weak effect of assumptions on IMFs. 



whose reported H magnitude implies a (lower limit) mass es- 
timate nearly one order of magnitude lower than J and K. We 
attribute this occurrence to an error of the integrated H magni- 
tude reported in 2MASS-6X as this value is at odds with that 
of all the other clusters while B347 is normal in all other re- 
spects. For instance it has a J-K color well within the range of 
the other clusters of the sample while its H-K color is more than 
one magnitude redder than any other. Finally we note that the 
independent lower limit mass obtained from the log(age) vs. My 



diagram (see Fig.[T4b. are in good agreement with that estimated 
from J and K magnitude for B347. Finally, as we have obtained 
just a lower limit to the age of B347 we do not provide an age 
estimate for this cluster. B347 as well as all the other clusters 
for which we can provide only a lower limit to the age are not 
included in the analysis of Sect. [5] that is limited to the young 
clusters that constitute the main subject of our study. 
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O MW OCs □ this survey 

X MW GCs |3 VDBO 




Fig. 14. Integrated V mag and total mass as a function of age for 
various samples of clusters. Galactic open clusters (OC, from the 
WEBDA database) are plotted as filled circles, Galactic globu- 
lar clusters (GC, My from the most recent version of the Harris 
(1996) catalog, i.e. that of February 2003, the ages have been ar- 
bitrarily assumed to be 12.0 Gyr for all the clusters) are plotted 
as x symbols. The target clusters are represented as open squares 
(VDBO as a crossed square), the clusters from P09b as open 
stars, and the clusters from WH01 clusters as open triangles. My 
magnitudes of the target clusters and of the P09b clusters are 
from the new aperture photometry performed on the CCD im- 
ages by Massey et al. (2006), except for B083 and B347 whose 
magnitudes are from RBC (see Tab. Q] My magnitudes of the 
WHOl's clusters are from RBC. Log Age is from Tab. [2] Points 
with arrows have only lower limits to the age. Filled circles are 
M31 OCs from Pap-II. The continuous lines are fixed-stellar- 
mass models from the set by Maraston (1998, 2005) for SSPs of 
solar metallicity, with a Salpeter's Initial Mass Function (IMF) 
and intermediate Horizontal Branch morphology. Note that in 
this plane, the dependence of the models from the assumed IMF, 
metallicity and HB morphology is quite small (see B08). The 
outlier OC at log Age^ 9.0 is Tombaugh 1. 



4. 1 . Comparison with Galactic open clusters 

In Fig. FBI we show the log(age) vs. absolute magnitude plot anal- 
ogous to Fig. Q~3]but using My instead of Mj, M H , M K . While 
NIR magnitudes are preferred to get reliable estimates of the 
stellar mass of our clusters (see Sect. [4] and Pap-I), the use of 
My allows us a direct comparison with different kinds of clusters 
for which integrated magnitudes in NIR passbands are lacking, 
Galactic OCs in particular (B08, Pap-I). 

Inspection of Fig. [14] confirms the tentative conclusions 
of Pap-I (and F05). The distribution of our target clusters 
marginally overlaps with the high-mass tail of the Galactic OC 
distributions, but the bulk of the sample of candidate YMC con- 
sidered here is significantly more massive than Galactic OCs 
in the same age range. In this sense, the brightest, most mas- 
sive and youngest cluster of our sample, VdBO having age=25 
Myr and M ^ 6 x 10 4 Mq, may appear similar to the handful 
of massive young clusters recently identified in the Milky Way 
(see Figer 120081 and Messineo et al. 120091 hereafter M09, for 
recent reviews), that have masses between 0.7 x 10 4 Mq and 
4.0 x 10 4 Mq and ages between 0.3 Myr and 18 Myr, according 
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X MW GCs £3 VDBO 

A WH01 




Fig. 15. Same as Fig. [T4]but with My magnitudes of the target 
clusters and of the WHOl's clusters obtained from fitting King 
(1966) models to our HST data, from Pap-III. The clusters from 
P09b are not included in the plot as they have not been consid- 
ered in Pap-III. 



to M09. The other clusters of our sample have similar (or slightly 
greater) masses than the Galactic YMC but they are all signifi- 
cantly older (by a factor of > 2x, see Sect. [5]for further discus- 
sion). It is worth to note that the masses estimated from Fig. [14] 
are in agreement with those from Fig.Q~3] typically, within a fac- 
tor of 2. 

In Pap-I we showed that in the case of VdBO, an exception- 
ally extended cluster, the integrated magnitudes reported in the 
RBC were significantly underestimated. However our shallow 
HST exposures were not ideal to perform integrated photome- 
try on such large areas (VdBO cover the whole extent of the PC 
field). For these reasons we recurred to the new homogeneous 
CCD survey by Massey et al. (2006; see Pap-I for discussion) to 
obtain a reliable estimate of the total luminosity of that cluster; 
as said, the integrated B,V magnitudes for the clusters consid- 
ered here have been obtained from the same source and with the 
same method (Tab Q]). These cases are less problematic, as the 
clusters are more compact than VdBO. However, it seems wise 
to check how the comparisons shown in Fig.[T4"lmav depend on 
the actual way in which My is estimated. To do that we present in 
Fig. [15] a new version of Fig.Q~4]in which the My values derived 
from Tab. Q] are replaced with My estimates obtained in Pap-III 
from profile fitting (with King 1966 models) performed on our 
HST images (with the same assumptions on distance and red- 
dening adopted here). Again, it is very reassuring to note that the 
conclusions drawn above from Fig. [14] are fully confirmed also 
by the new set of My from Pap-III. In fact, the differences be- 
tween the YMC of our sample and Galactic OCs are even more 
pronounced in the new plot, as the total V luminosities estimated 
in Pap-III are larger than the values adopted here by a factor of 
^ 1.6, in average. For the reasons discussed in Pap-I and for 
homogeneity with that analysis we retain our ground-based My 
estimates as our reference. 

It is interesting to note that the clusters identified by Krienke 
& Hodge ([2007112008b , and, by analogy, those found in Pap-lfl 



7 It should be recalled that clusters listed in the RBC were excluded 
from the analysis performed in Pap-II. 



S. Perina et al.: A HST/WFPC2 survey of bright young clusters in M31. II. 



19 



20 

: 22 
24 



-■ 1 1 1 1 1 


1 1 1 


-' 1 1 1 1 1 


1 1 1 '- 


- t=251Myr 


B448 - 


- t=316Myr 


B081 - 


- Z=0.03 




- Z=0.03 




E(B-V)=0.55 




E(B-V)=0.28 












-, I , 1 , I 


, , 1 


" , 1 1 , , 1 


, , 1 



2 4 

F450W-F814W 



2 4 

F450W-F814W 



is 

"3 
u 



T 8.5 



< 
o 



I I I I I I I I I I I I I I I I I I ..I I I /I 




7.5 8 8.5 9 9.5 
log Age [yr] (this survey) 

Fig. 16. Bottom panel: comparison of the CMD-based ages from 
Tab. 2 with the ages obtained by C09 from integrated spectra. 
The symbols are the same as in Fig. [14] B257D is not plotted be- 
cause it is not included in the C09 sample. The error bars show 
the average errors. The vertical arrows indicate clusters defined 
as "older" than 2 Gyrby Caldwell et al. (2009). The two clusters 
from our own survey for which the two independent estimates 
show the greatest difference are labeled (B448 and B081). Top 
panels: Comparison of the observed CMD for B448 and B081 
with the isochrone corresponding to the age, metallicity and red- 
dening estimates provided by C09 for these clusters (values re- 
ported in the upper left corner of each panel). Note that in the 
case of B448 the reddening estimated by C09 is obviously too 
large, while in the case of B08 1 , the metallicity assumed by C09 
(Z=0.03 for all the clusters) seems the principal responsible for 
the mismatch. 



have an observed LF peaking around My = -3 and virtually 
dropping to zero at My ^ -6, very similar to Galactic OCs (see 
Fig-QSK hence they appear as the natural counterpart of the OCs 
observed in the Milky Way. 

In Pap-III the problem of the survival of our target clusters 
was discussed in some detail and dissolution times including 
the effects of internal and external evolution (Lamers & Gieles 
2006), were computed. These values are reported also here, in 
Tab. [3] for convenience of the reader. The dissolution times of 
young clusters are all shorter than a Hubble time, hence it is 
likely that none of them will survive long enough to become 
old (age^ 10 Gyr), and some of them are probably in the lat- 
est phase of their dissolution (B321, B342; Pap-III). However, a 
few clusters have dissolution times longer than 1 Gyr, and it is 
not inconceivable that some of them may reach an age of several 
Gyr before dissolving into the M31 disk (see Pap-III). 



4.2. Comparisons with Caldwell et al. (2009) 

A comparison of the results obtained here from the analysis of 
our HST-WFPC2 CMDs with those of the extensive and the 
independent analysis by C09, based on high-quality integrated 
spectra is clearly worthwhile, in this context. 

In the lower panel of Fig. [16] the age estimates from Table. [2] 
are compared with those by C09. The two set of ages do agree 
within the uncertainties, but there is a clear systematic offset as 
C09 ages are larger than those listed in Tab. [2] by a factor of 
s 1.5, in average, and up to a factor of ^ 3 in the worst case 
(we are considering only clusters having age estimates in both 
sets, not lower limits). We note that this systematic offset occurs 
also if one restricts the sample by WH01, and also to the three 
clusters for which C09 provides CMD-based age estimates of 
their own (see their Tab. 7), hence it is a characteristic feature of 
their spectroscopic age estimates. 

A difference that may produce a systematic offset between 
our ages and those by C09 is that they adopt super-solar metallic- 
ity models (Z = 0.04) for all the clusters, while we leave metal- 
licity as a free parameter of our fit and, in fact, we adopt solar 
or less-than-solar metallicity models in all cases (see Tab.FJ}. If 
both sets of ages were derived from isochrones fitting the effect 
should be the opposite, i.e. a younger isochrone is required to 
fit a given CMD with a model of higher metallicity. However it 
is not clear if this general behavior is shared also by models of 
integrated spectra. 

In the upper panels of Fig.[16]we show the two cases (among 
those included in our own survey) that display the widest dif- 
ference between the two age estimates. We superposed on the 
observed CMDs the isochrones corresponding to the best-fit es- 
timates by C09, corrected by the reddening provided by these 
authors. The case of B448 shows very clearly that the solution 
provided by C09 significantly overestimates the reddening, and 
it is not compatible with the observed CMD. In the case of B081, 
the comparison suggests that the choice of super-solar metallic- 
ity models by C09 may be particularly unsuitable for this cluster, 
leading to a larger-than-average error in the age estimate. 

Two cases of especially remarkable differences occur also 
with the set by WH01 (open triangles in Fig. [16}. B319=G44 
is considered also in Tab. 7 of C09, where a spectroscopic age 
of 0.28 Gyr is reported, to be compared to the CMD-based age 
estimated of 0.10 Gyr by WH01. Moreover the reported spec- 
troscopic value is most probably a typo, as in Table 2 of C09 
(their primary source of cluster ages) they report log(age)=8.6 
for B319=G44, corresponding to 0.398 Gyr (the value that is 
plotted in Fig. [16]). In any case, the spectrum appears to be rea- 
sonably fitted by a Z=0.04, age=500 Myr model (N. Caldwell, 
private communication), while the CMD shown by WH01 is 
clearly not compatible with such an old age. The a-priori as- 
sumption of super-solar metallicity models by C09 may also be 
the origin of this mismatch. The case of B368=G293 (not in- 
cluded in Tab. 7 of C09), that is classified by C09 as "older than 
2 Gyr" while the CMD by WH01 indicates age ^ 80 Myr, has 
to be ascribed to a typographical error by C09; in fact the cluster 
was not observed by that authors (N. Caldwell, private commu- 
nication). 

Fig.fTTIshows the comparison between our estimates of E(B- 
V) and those by C09. In this case as well there is reasonable 
overall agreement, most of the differences being within the un- 
certainties. The most discrepant case is B448, already discussed 
above (see Fig. [TBI). Finally, in Fig. [18] the mass estimates are 
compared. Also in these cases the two set of estimates agree 
within the uncertainties (1 cr is a factor of 2.4), the strongest 
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Table 3. Newly derived masses and dissolution times for the studied clusters. 
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In a few cases the data allowed us to obtain only a lower limit to the mass, tt/2000 and djiooo are from 2MASS-6X-PSC catalog, J, H, K are from 
r=1070 ap. phot, in the 2MASS-6X-PSC catalog. Note that err JHK = l.OO corresponds to err JHK =nu\l in the 2MASS-6X-PSC catalog. 
Estimated from Fig. [14] as these clusters lack NIR photometry. These mass estimates will not be used in the following to preserve the 
homogeneity of the sample. 



discrepancy is to be attributed to the overestimate of the age for 5. Summary and Discussion 

B319=G44 by C09 discussed above. 

We presented the main results of a survey aimed at the determi- 
nation of the nature of a sample of 20 candidate YMC in the thin 
disk of M31 (one of which, VdBO, was studied in Pap-I). One of 
the targets surveyed turned out to be a bright star projected onto 
In conclusion, while we are unable to identify the reason of the dense disk of M31, and thus erroneously classified as a pos- 
the (modest) systematic overestimate of the ages by C09, it has sible cluster. All the other targets were revealed to be genuine 
to be concluded that the agreement between the two independent star clusters and we were able to obtain reliable CMDs for all of 
sets of age, reddening, and mass estimates is quite satisfactory, them. The main results from our own survey can be summarized 
if the observational uncertainties are taken into the due account, as follows: 
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Fig. 17. Comparison of the E(B-V) estimates from Tab. 3 with 
those by C09. The symbols are the same as in Fig. [14] 
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Fig. 18. Comparison of the masses estimates from Tab. 3 with 
those by C09. The symbols are the same as in Fig. [14] The grey 
symbols show the clusters that have "null" error on IR mag- 
nitudes in the 2MASS-6X-PSC catalog. The thick line is the 
Mt. s . = Mcq9 locus, the thin lines bracket the + l<x range about 
this locus. The error bars show the average errors. 



1. New integrated-light spectroscopy became available for 
many of our targets since the original selection was per- 
formed. Three of them (B083, NB16 and B347) were re- 
vealed by the new data to be not good YMC candidates as 
defined by F05. The CMDs obtained in this study confirms 
that they are likely old clusters. 

2. Among the remaining 17 targets, 16 are genuine clusters 
and one is in fact a star (NB67), as said above. Thus the 
fraction of spurious objects in our well-defined sample of 
BLCC=YMC is just 1/16 = 6.2%. Even excluding the two 
clusters considered at point 3., below, the incidence remains 
below 10%. The extended sample considered in AppendixlBl 
fully confirms these results. We must conclude that M31 
YMC are not especially plagued by contamination from spu- 
rious sources and most of the clusters considered in the orig- 
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Fig. 19. Upper panel: The mass distribution of the sample of 
YMC studied here (from Tab. [3] thick continuous line) is com- 
pared with the mass distribution of Galactic OCs (dotted line) 
and Galactic globular clusters (dashed lines). Masses of Galactic 
clusters are from B08. Lower panel: zoomed view of the distri- 
bution of M3 1 YMC compared with the distribution of the YMC 
of the Milky Way (data from M09). 

inal analysis by F05 should be rea0- In particular, asterisms, 
suggested as a possible major contaminant of the sample by 
C06, are in fact found to be not a particular reason of con- 
cern, in this context (see also the discussion by C09). 

3. Two of the sixteen genuine clusters (B374 and B222) have 
integrated properties compatible with being YMCs but they 
do not show a detectable MS in the range of magnitudes sam- 
pled by our CMDs. We can provide only an upper limit to the 
age of these clusters (^ 300 Myr), but the available data sug- 
gest that they are good candidate intermediate-age clusters 
that indeed would merit follow-up with deeper HST photom- 
etry. 

4. The fourteen confirmed young clusters (including VdBO, 
studied in Pap-I) show a clear MS in the range of magni- 
tudes sampled by our CMDs, hence we were able to obtain 
reliable estimates of their ages, reddenings and (an educated 
guess of) metallicities by comparison of the observed CMD 
and LF with theoretical models. Ten of them have ages in the 
range 25 - 1 00 Myr, the other four range between 1 40 Myr and 
280 Myr. The adopted metallicities include Z = 0.004 (one 
case), Z = 0.008 (three cases), and Z = 0.019 (solar metal - 
licity, ten cases). The estimated reddenings range from E(B- 
V)=0.06 to E(B-V)=0.60, with E(B-V)=0.20-0.30 as most 
typical values. 

To increment our final sample of YMC we included ten fur- 
ther clusters for which the age was estimated from their CMDs 
(obtained from HST imaging) with methods strictly homoge- 



8 It may be useful to stress again that the clusters of our survey were 
selected among the class f=l RBC entries, see Sect. [2] and Galleti et 
al. d2006at . 
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Fig. 20. Comparison between Galactic OCs (small filled circles), M3 1 YMC from the present study (big open squares), MW YMC 
from M09 (big open circles), M3 1 's clusters from pap-II (small open squares), Magellanic Clouds clusters (grey open pentagons), 
and M33's clusters (grey crosses) in the log(age) vs. log Mass plane. Masses of Galactic OCs are from B08, masses of Magellanic 
Clouds clusters are from 2006 and masses of M33 clusters are from 2009, For M33 and the Magellanic Clouds only clusters younger 
then 10 Gyr are shown. 



neous with those adopted here, from WH01 and P09b. In this 
way we assembled a final sample of 24 confirmed young clus- 
ters. For 22 of these we were able to obtain reliable estimates 
of the total stellar mass by coupling our age estimates with the 
integrated J,H,K magnitudes taken from the 2MASS-6X cata- 
log. These clusters have masses ranging from 0.6 x 10 4 M© to 
6 x 10 4 M Q , with an average of ~ 3x 10 4 MeH Our estimates of 
ages and masses are in good agreement with recent independent 
studies based on integrated light spectra (see also Pap-III for the 
comparison with the results by Pfalzner 20d5)l. 



9 The remaining two clusters, that lack NIR photometry, also have 
masses lying in the same range, according to the estimates obtained 
using the integrated V magnitude instead of J,H,K ones. 



5.1. The nature of M31 YMC 



In the upper panel of Fig. [19] the mass distribution of our ex- 
tended sample of M3 1 YMCs is compared with the distributions 
of Galactic OCs and GCs (masses from B08). The clusters con- 
sidered here appear to lie in the middle of the two distributions, 
overlapping with the high-mass end of the OCs and with the low- 
mass end of GCs. This comparison provide a further confirma- 
tion that the YMCs (=BLCCs) of M31 are indeed more similar 
to the YMCs of the LMC than to classical OCs of the Milky 
Way, i.e. the original hypothesis advanced in F05. This is in full 
agreement with the main conclusions by C09, obtained with a 
completely independent method (less sensitive to age than ours) 
on a wider sample. 
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The lower panel of Fig. [19] compares our clusters with the 
YMCs seen toward the center of the Milky Way as listed by M09. 
The two samples have very similar mass distributions, suggest- 
ing that they are also similar in nature. An obvious difference 
between the two sets of clusters was already suggested in Pap-I 
and is confirmed here: the M31 YMCs of our sample are signif- 
icantly older that the YMC discovered until now in the Galaxy 
(k. 50 Myr vs. $ 20 Myr; see below for possible explanations). 
We confirm that the M3 1 YMCs studied here have larger sizes 
(half-light-radii) with respect to their MW counterparts (see Pap- 
I and Pap-III); this seems in agreement with the age-size rela- 
tions proposed by Pfalzner (2009; see Pap-III for discussion). 

A more thorough comparison between various samples of 
YMCs is presented in Fig. [20] where Galactic OCs and YMCs, 
YMCs from M33 (San Roman et al. [20091 for further discus- 
sion on M33's star clusters see Sarajedini & Mancone 2007, 
Zloczewski et al. [20081 Park et al. |2U09l , the LM C the Small 
Magellanic Cloud (McLaughlin & van der Marel[2006), and 
M31 are plotted together in a log(age) vs. log Mass diagram. 
Fig. [20] is affected by a number of selection effects that deserve 
to be described in some detail. 

1 . The minimum mass threshold appears to increase with age 
(at least for age ^ 10 Myr, see the Galactic OCs if Fig. [20b: 
this is due to the fact that the lower the mass of a cluster, the 
shorter is its dissolution time, as the cluster is less resilient 
to all the internal and external effects that may lead to its dis- 
ruption (Gieles et al. 120071 Pap-III, and references therein). 
The minimum mass threshold for samples in external galax- 
ies is obviously due to the inherent magnitude limits. 

2. Also the maximum mass threshold increases with age in 
log Age vs. log Mass plots (Hunter et al. 120031 Gieles |2009l 
the effect is clearly evident in Fig.[20]if one looks at the MW 
OCs, that cover the widest range in ages). This general be- 
havior can be easily explained as a simple consequence of 
varying the sample size as a function of the age bin in the 
logarithmic scale. Assuming a power-law mass function and 
a constant Cluster Formation Rate (CFR) the number of clus- 
ter per logarithmic age bin increases with age. For an expo- 
nent of the power law mass function (N(M) oc M~ a ) a — 2, 
that is a reasonable approximation for most of the observed 
cluster systems, log M max oc log Age (see Gieles 2009] for 
detailed discussion and references). 

3. While the lack of massive (M ^ 10 4 Mq) clusters older than 
400 Myr in the Milky Way is probably real, the typical limit- 
ing magnitude (V ~ 27, Rich et al. 2005) of available CMDs 
of M3 1 clusters prevent us from drawing firm general con- 
clusions about objects in that age range in M3 1 . The cases of 
B222 and B374, treated here, are excellent examples of clus- 
ters that may populate that region of the diagram but lack a 
reliable age estimate because the available photometry is too 
shallow (see Puzia et al. |2005l ). 

4. The lack of massive (log (M/Mq) > 3.6) M31 clusters 
younger than 25-50 Myr may be due to the contribution 
of several biases. First, such young clusters may be hard 
to select from the RBC as there are no objects bluer than 
{B - V)o - 0.0 in the list of confirmed clusters (see F05). 
This is not surprising as the RBC was intended to be a cat- 
alog of globular clusters. Second, for ages SS 8 Myr the Hp 
index is expected to fall below the threshold adopted to select 
YMC candidates (see, for example, Fig. 7 of F05), thus (pos- 
sibly) preventing the selection of these objects for our survey. 
Third, very young objects should have their luminosity dom- 
inated by a few massive stars near their centers, thus leading 



to objects that may appear more like blended stars than like 
a star cluster at the distance of M31, even in HST images, 
thus preventing their inclusions in lists of candidate YMCs. 
Fourth, it can be hypothesized a positive correlation between 
the age of the clusters and their height above the disk plane, 
such that the youngest clusters are more deeply embedded 
in the thin dust layer of the M3 1 disc, out of our reach even 
from our privileged point of view, while most/some of the 
older clusters would be visible just because they lie above 
the densest part of that layer. There are indications that this 
kind of correlation actually holds in our own Galaxy (V.D. 
Ivanov, private communication). 

5. The lack of massive (log (M/Mq) > 3.6) MW clusters older 
than 25-50 Myr may also be associated with an observational 
bias. Galactic YMC have been identified as clumps of bright 
stars in the near and mid IR and the youngest clusters, having 
the brightest RSG, are easier to detect in this way. Moreover 
the sample of Open/YM Galactic clusters is limited (essen- 
tially by the effect of interstellar extinction in the Galactic 
disc) to a volume of a few kpc around the Sun, while M3 1 
(or M33) YMCs can be selected over the whole disk of their 
parent galaxy, thus introducing a bias that favors the detec- 
tion of rarer cluster species (massive clusters) in the latter 
galaxies with respect to the MW. 

6. There seems to be a significantly under-dense region in 
Fig. 1201 for masses <: 10 3 Mq and ages between ~ 15 Myr 
and ~ 50 Myr (7.2 i= log Age^ 7.7). The same feature was 
noted by Whitmore et al. (120071 1 in their study of the clus- 
ter system of the Antennae and it was attributed by a de- 
generacy in age dating from broad band colors occurring in 
that age range due to the prompt onset of the RSG phase 
(see Whitmore et al. 2007 , for details, discussion and further 
references). Virtually all the clusters plotted in Fig. [20] had 
their ages estimated from the CMD of their stars (instead of 
broad-band colors, see also Pap-II), hence our sample should 
not be affected by this bias, at least in principle. However the 
coincidence of the feature with that noted by Whitmore et 
al. (12007) suggests that the same kind of bias against ages in 
that interval may be at work also in Fig. [20] 

7. The samples of clusters from all the galaxies involved in 
Fig.[20]have been selected according to different criteria, by 
color, magnitude, etc. 



Given all the above considerations, it does not seem possi- 
ble to draw any firm conclusion from the comparison shown in 
Fig. [20] The only straightforward conclusion is that YMCs in 
the age range 50-500 Myr are relatively common in all the most 
massive star-forming galaxies of the Local Group (M31, M33, 
LMC and SMC). The only exception (the Milky Way) may be 
ascribable to observational biases, but it cannot be excluded that 
it is instead (at least partly) associated with intrinsic properties of 
the Milky Way, that appears peculiar under several aspects with 
respect to the typical spiral galaxies (and to M31, in particular 
see Hammer et al. 2007 and Yin et al. |20091 >. As the samples 
of M33 and M31 should be subject to the same kind of biases 
(as the distances are similar and the data have been collected 
with HST in both cases), the difference in the maximum mass 
limit between the two samples is likely real, and it can probably 
be ascribed to the difference in total mass between the discs of 
the two galaxies: larger discs should host more numerous popu- 
lations of clusters, thus enhancing the probability of producing 
clusters with higher (maximum) masses (see Gieles 2009, and 
references therein). 
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5.2. Radial trends 

Given the wealth of data collected for our target clusters, it may 
be useful to look for correlations between their physical param- 
eters, including their position within the M31 disc. Limiting the 
analysis to the young clusters (age < 1 Gyr), that constitute 
a more homogeneous sample of bona-fide thin disk objects, it 
turns out that our sample is still too sparse for a thorough analy- 
sis of these correlations. In particular the covered ranges of age, 
mass and position are quite limited, thus not allowing us to re- 
veal large scale trends, in most cases. Moreover, the adopted ap- 
proach of CMD analysis provides just an educated guess of the 
metallicity of the clusters, aimed at obtaining the most reliable 
estimate of the clusters age, which was the main objective of 
our analysis. These limitations prevent the possibility of a mean- 
ingful study of the radial metallicity gradient with our data. It 
should also be recalled that the correlations bewteen the struc- 
tural parameters of the clusters (mass, radius, density etc.) have 
already been discussed in Pap-III, hence here we consider only 
age, mass, de-projected galactocentric distance (Rj\ assuming 
and inclination of i = 12.5° of the disk with respect to the plane 
of the sky, see Simien et al. 1978 and Pritchet & van den Bergh 
[19941 , X, Y, and reddening. 

Having checked all the combination of parameters, the 
only correlation that appeared remarkable to us is presented in 
Fig. [2T] It is a trend of decreasing age with galactocentric dis- 
tance, that seems statistically significant if one consider the as- 
sociated errors. Given the relatively limited range of galactocen- 
tric distance covered, in our view the observed distribution can 
be interpreted in two ways: 

- as a part of a larger trend resulting from a inside-out wave of 
cluster formation. In this case the trend toward older mean 
ages should continue at lower radii and Fig. [21] shows the 
transition between a regime of decreasing age with galacto- 
centric distance and an asymptotic regime of constant age in 
the outermost fringes of the disc; 

- more likely, as a sharp transition in the epoch of the highest 
rate of star/cluster formation occurring at the onset of the 
Rd ~ 10 kpc "ring of fire". This would be consistent with the 
well known burst of recent star formation that characterize 
this prominent structure of the M3 1 disc. 

While not especially conclusive or insightful, the result 
shown in Fig. [21] gives a clear idea of how useful YMCs can 
be as tracers of the structure and evolution of the disk itself, in 
particular if large and reliable samples can be assembled. 

5.3. Final remarks 

This research has demonstrated that the conspicuous popula- 
tion of bright disk objects studied by F05 consists of genuine 
YMC, similar to those found in the LMC, SMC and M33 galax- 
ies. These clusters may open a new window to the study of the 
recent star formation history in the disk of M3 1 . A systematic 
analysis over the whole extent of the M31 disk may provide 
the opportunity to study a rich system of young clusters using 
a sample much less affected by selection biases than in our own 
Galaxy, and to better constrain the models of dynamical evolu- 
tion of clusters within the discs of spiral galaxies. M31 YMCs 
like those studied here provide also an excellent tracer of the disk 
kinematics in that galaxy, independent of (and in addition to) the 
HI gas. Recent wide-field surveys (Vansevicius et al. 20091 see 
also Pap-II) suggest that a rich harvest of genuine YMCs await 
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Fig. 21. Age as a function of the deprojected galactocentric dis- 
tance for the young clusters (open squares with error bars). The 
cluster VdBO has been labeled as it is by far, the youngest of the 
whole sample. 

to be discovered in the disk of our next neighbor giant galaxy in 
Andromeda. 
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Appendix A: RBC clusters serendipitously imaged 
in our survey 

To ascertain the real nature of candidate M3 1 clusters proposed 
by various authors is a daunting but necessary task to keep clus- 
ter catalogs as complete and clean as possible from spurious 
sources. There are several criteria that may be used to check can- 
didates (see Galleti et al. 2006a for references and discussion), 
but resolving them into stars by means of high spatial resolu- 
tion imaging is by far the safest method of all. In addition to the 
clusters that were the main target of our survey, and to the low- 
luminosity clusters identified by Hodge et al. 2009} our WFPC2 
images serendipitously included several clusters and candidate 
clusters listed in the RBC. Inspection of our images allowed us 
to place their classification on firmer footing. The results of this 
analysis are summarized in Table I A. fl Their classification in the 
RBC has been modified accordingly. In Table lA.fl we report the 
name of the object (column 1, name), the classification flag orig- 
inally reported in the RBC (col. 2, f), the name of the cluster that 
was the original target of the images (col. 3, field), a flag indicat- 
ing if the object was imaged with the PC or with one of the WF 
cameras (col. 4, chip), and, finally, a comment on its classifica- 
tion as derived from the inspection of the new images. In some 
case the classification remains uncertain (comments with "?"). 
In some cases the image reveals that the object is extended but 
do not clarify its nature (cluster/galaxy/HII region etc.), in these 
cases we report the comment "not a star". An estimate of the ra- 
dial velocity will suffice to definitely establish if these objects are 
M31 clusters or background galaxies (see Galleti et al. [2006a). 
In some cases, some clusters that were among the main targets 
of our survey were serendipitously re-imaged in the WF field 
surrounding other targets. For obvious reasons these cases are 
not reported in Table I A. f| On the other hand some clusters have 
been serendipitously imaged in two different pointings: in these 
cases we report the classification derived from both sets of im- 
ages. Some of the clusters of Table |A.1~| were independently re- 
identified in Pap-II (B061D, B319, B014D, B256D, DA084), for 
two of them a meaningful CMD was also obtained there (B061D 
and B319); this lends additional support to the reliability of their 
classification. Finally, we reported in the table also some clus- 
ters whose nature was already confirmed by previous HST imag- 
ing, for completeness (see the case of B3 19=G044, observed by 
WH01). 

It may be interesting to note that among the 19 RBC class 
f=2 (candidate clusters) objects listed in Tab. IA.1I 3 turn out to 
be real clusters (or likely clusters), 5 are extended objects that 
lack the v r measure needed to ultimately establish their member- 
ship to M31, while 11 are non-clusters (or likely non-clusters), 
most of them being stars. According to this limited sample it can 
be concluded that the fraction of genuine M3 1 clusters among 
class f=2 entries of the RBC ranges from -j| = 16%± 14% to 
y|=42%+ 12%. These numbers should be considered as some- 
what pessimistic as they are computed on a sample of clusters 
projected on the densest regions of the M31 disc, where the 
probability of contamination from bright stars of M3 1 is at its 
maximum. To give a rough idea of the number of genuine clus- 
ters that are still hidden among the candidates listed in the RBC 
one can take the 16% of the number of class=2 RBC entries, i.e. 
0.16 x 1049 168. A significant fraction of these may be YMCs 
(> 15%, according to F05). 

Considering the objects listed in Tab.Q]and Tab. [AT] the sur- 
vey images allowed us to verify the nature of 25 objects classi- 
fied as genuine clusters (class f=l) in the RBC. We confirm that 
23 of them are real clusters while 2 are (one or two) stars. From 



this number one can estimate the fraction of spurious sources 
among class f=l RBC entries as ^=8%+8%, that is remarkably 
low and is in excellent agreement with the estimate by G09 that 
finds 5=4% from a sample of 252 objects. 

Considering the fraction of real clusters among class f=l en- 
tries as 92% and that among f=2 entries as 16%, the expected 
number of genuine M3 1 clusters in the RBC (GC+ YMC) is es- 
timated as ~ 630, while the number of old clusters (GCs) should 
be ~ 530, in reasonable agreement with the results by Barmby et 
al. 2000 and F05. Note that, at present, the number of confirmed 
(likely) old clusters (f=l and y=0) in the RBC is 418; correcting 
this for contamination leads to 384 bona-fide GCs, more than 
double than the number of GCs encountered in the Milky Way 
galaxy (=* 150, Harris [19961 

Appendix B: Other candidate M31 YMCs with 
archival HST imaging 

Before selecting the actual targets for our survey we searched 
the HST archive for YMC candidates, as listed in Tab. 1 (or 
Tab. 2) of F05, that had already been (serendipitously) imaged 
from HST. As the nature of these objects (cluster / asterism / 
star) can be determined from existing images they were not in- 
cluded in our final list of targets. In Tab. B.l. (referring to ob- 
jectively selected candidates from Tab. 1 of F05) and Tab. B.2. 
(referring to candidates suggested from various authors adopting 
different criteria, from Tab. 2 of F05) we list the results of that 
research. In these tables we report (1) the cluster name(s), (2) the 
HST program number(s) of the retrieved images, (3) the instru- 
ments) and (4) the filter(s) used to obtain the inspected images, 
(5) the classification of the object based on the inspection of the 
HST images, following the approach adopted in Tab. IA.ll above, 
and, finally, (6) the classification provided by C09 based on their 
spectra and/or on ground-based imaging (S indicates that the ob- 
jects was classified by from its spectrum, I indicates that the ob- 
ject was classified with imaging, SI means that both imaging 
and spectrum were considered for the classification, according 
to C09). At the epoch when the table was compiled (September 
2009), 36 out of the 66 objects listed in Tab. 1 of F05 (including 
those studied in this paper) had one (or more) images in the HST 
archive: 34 of them are recognized as real star clusters from the 
inspection of the available HST images, while 2 are stars. This 
leads to a fraction of spurious objects in the sample of 5.5% 
+ 4.0%, in full agreement with the fraction we obtained from 
our original sample (Sect. [2]). Analogously, 14 out of 21 objects 
listed in Tab. 2 of F05 (including those studied in this paper) had 
one (or more) image(s) in the HST archive: 13 of them are rec- 
ognized as real star clusters from the inspection of the available 
HST images, while 1 is a star. This leads to a fraction of spurious 
objects in the sample of 7.1% + 7.4%, again in full agreement 
with the fraction we obtained from our original sample (Sect. [2]) 
and with the above results. Note that (a) all the classifications 
we obtained from HST imaging confirm those independently ob- 
tained by C09 for the same objects, and (b) all the objects listed 
in Tab. B.2. were classified as clusters by some other author be- 
fore (see F05). 

Of the 37 objects in Tab. B.l. and Tab. B.2. lacking HST- 
based classification, 31 are classified as clusters by C09; the 
remaining 6 have uncertain classification. Coupling the results 
from HST and C09 it turns out that 60 of the 66 objects from 
Tab. 1 of F05 are real clusters, two are stars, and four have un- 
certain classification; 18 of the 21 objects from Tab. 2 of F05 are 
real clusters, one is a star, and two have uncertain classification. 
We thus conclude that the large majority (^90%) of the objects 
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identified (or proposed) by F05 as (possibly) young clusters are 
indeed genuine star clusters. Finally, three clusters listed in the 
RBC but not comprised in the study by F05 where found in Pap- 
II to have age < 1 Gyr (B014D, B061D, B256D). 
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Table A.l. RBC clusters serendipitously imaged in our survey. 



Name 


f 1 


Field 


Chip 


Comment 


B014D 


2 


B015D 


PC 


cluster 


B061D 


2 


NB16 


WF 


cluster 


B256D 


2 


B257D 


WF 


cluster 2 


B256D 


2 


B475 


WF 


cluster 2 


SK067B 


2 


B015D 


WF 


not a star 


SK071C 


2 


B475 


WF 


not a star 


SK185B 


2 


B475 


WF 


not a star 


B068D 


2 


NB16 


WF 


not a star 


B068D 


2 


NB67 


WF 


not a star 


B019D 


2 


V031 


WF 


not a star 


NB64 


2 


NB16 


WF 


star? 


NB64 


2 


NB67 


WF 


star? 


SK091B 


2 


B066 


WF 


star 


B048D 


2 


B081 


PC 


star 


SK091C 


2 


B374 


WF 


star 


SK188B 


2 


B475 


WF 


star 


NB47 


2 


NB16 


WF 


star 


SK083B 


2 


B043 


WF 


2 stars + nebula? 


B057D 


2 


NB16 


WF 


2 stars 


NB43 


2 


NB67 


WF 


2 stars 


B192D 


2 


B327 


WF 


galaxy 


SK194C 


2 


B376 


WF 


galaxy 


B376 


1 


B374 


WF 


cluster 


B257D 


1 


B475 


WF 


cluster 


B319 


1 


B318 


WF 


cluster 


DA084 


1 


B374 


WF 


not a star 3 


DA084 


1 


B376 


WF 


not a star 3 


SK047A 


1 


B081 


WF 


two stars 


NB68 


6 


NB16 


WF 


star? 


NB68 


6 


NB67 


WF 


star? 


B113 


6 


NB16 


WF 


star? 


SK069D 


6 


B083 


WF 


star 


B185D 


6 


B318 


PC 


star 


SK046D 


6 


B327 


WF 


star 


B065D 


6 


NB67 


WF 


star 


SK041D 


6 


B321 


WF 


two stars 


B121 


3 


NB16 


WF 


star? 


B121 


3 


NB67 


WF 


star 



1 f is the original RBC classification flag (1 globular cluster, 2 candidate globular cluster, 3 controversial object, 6 star/s). 

2 While the visual inspection of the images does not permit a clear cut classification, the objective analysis performed in Pap-II recognizes 
B256D as a star cluster. 

3 DA084 has a radial velocity estimate that clearly identifies it as a member of M31 (see the RBC). 
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Table B.l. Classification of candidate young clusters listed in Tab. 1 of F05. 



Name 


Obs-ID 


Camera 


Filters 


Class HST 


Class C09 


BUUS-GUoU 


10407 


Ape Ai/cr 1 


F606W F43jW 


cluster 


cluster(SI) 


BUZfi-GUoS 










cluster(SI) 


DU4U-L1UZ 


108 18 


WrrLZ 


C/1 «A\li/ CO 1 AW! 

r4jUW roi4W 


cluster 


cluster^!) 


B043-G106 


10818 


WFPC2 


P/1 CAM/ POl AW7 

F450W rsl4W 


cluster 


cluster(SI) 


D A/l"7 i" 1 1 1 1 

BU4/-G1 1 1 










cluster(S) 


D A /I A i" 1 1 1 O 

B049-G1 12 


1040 /( 10631) 


ALcy WrL 


F435W ro06W 


cluster 


cluster (SI) 


B057-G1 18 


10407(10631) 


ACS/WFC 


F435W F606W 


cluster 


cluster (SI) 


BU66-G1Z8 








cluster 


cluster(M) 


D Ait A 1 -JO 
DUby-Lll JZ 


10273 


AL5/WrL 


PjjjW rol4W 


cluster 


cluster(SI) 


B074-G135 










cluster(S) 


D AO 1 1~> 1 /lO 


10818 


WrrLZ 


C/1 Zt\WI CO 1 AWT 

r4jUW rol4W 


cluster 


cluster(ol) 


D A 1 i" 1 1 /I <i 

BU83-G146 


10818 


WrFLz 


C/1CA"\17 CO 1 AWT 

r45UW r8l4W 


cluster 


cluster(S) 


D AA 1 PK1 

BUyi-Gl M 


10273 


APC Al/up 

ALi/WrL 


PjjjW rol4W 


cluster 


cluster (SI) 


D 1 14-U1 / J 


5907 


WrrLZ 


C<^ <WT CO 1 AWT 


cluster 


cJuster(ol) 


rS 16U-UZ14 


fl/IOA/l A1T3 "7 A^C^\ 

y4oU( luz / d, /4Zo ) 


ALo/WrL, WrrLZ 


C"7"7 ZWT CZ^ZWT CO 1 AWT C£A£Y!l' 

r / / j W rOjj W rol4 W roUoW 


cluster 


cluster(ol) 


d 1 /U-UzZl 










.-.1,, C FA 

cluster(M) 


dZIU-MI 1 


9709 


WrrLZ 


F606W 


cluster 


cluster(SI) 


BZ16-GZ6/ 










cluster(SI) 


dZZZ-KjZ 1 1 


10818 


WrrLZ 


C/1 5n\X7 CO 1 AWT 

r4jUW rol4W 


cluster 


cluster(ol) 


B223-G278 










cluster (SI) 


Din i~">aa 

bzj /-Gzyy 










cluster(Sl) 


rSZol-LZoo 










cluster (SI) 


Dine i" 1 A 1 /I 

Bzy5-GU14 










cluster(S) 


D1AQ r^CiOA 










cluster(SI) 


B30/-GU30 










cluster(SI) 


B314-GU3 / 










cluster(^l) 


nil c o 
B315-GU3fi 


8296 


WrrLZ 


C11^\1! C/11A\i/ CCCCA17" 

F336W F439W F555W 


cluster 


cluster (SI) 


DUO ("'A/IO 
djI o-LU4Z 


oZy6(lufilfi) 


WrrLZ 


cii£.wt c/i invir tja ^ a\s/ kc^ wt co 1 a wi 
rjJ6W r4JyW r4jUW rjjjW rol4W 


cluster 


cluster(^l) 


duo t~*f\A A 
dj iy-LjU44 


8296 


WrrLZ 


Cllt.Wf CAKIWT CA^I\WT C^^^WT CQ1/1W/ 

F336W r43yW r4jUW FjjjW rfi!4W 


cluster 


cluster(SI) 


B321-GU46 


10818 


WFPC2 


P/1 CAM/ rOI AWT 

F45UW F814W 


cluster 


cluster(SI) 


Din pAlA 

B3zz-Gu4y 










cluster(^l) 


B32/-GU53 


10818 


WrrLZ 


C/1CA\SF CO 1 AWT 

F450W r8l4W 


cluster 


cluster(SI) 


B331-G057 


6699 


WFPC2 


F555W F814W 


cluster 


cluster (SI) 


B34z-GUy4 


8296 


WrrLZ 


ctj^uj c/nrnu rcccui 

F336W F43yW F555W 


cluster 


cluster(^l) 


B3j4-G1<S6 










cluster(S) 


B355 










possible star(S) 


D1CQ f) 1 A 

d j jo-oz iy 










candidate 


B3o/-Gzyz 


10407 


ALb/WrL 


F435W roOoW 


cluster 


cluster(SI) 


t> jofi-ozyj 


ozyo 


WrrLZ 


cn^w c^^w 
rjjOW r4jyw rjjjW 


cluster 


cluster(I) 


D1"7/1 PQAA 

dj /4-OjUD 


108 18 


WrrLZ 


C/1 ^A\S/ CO 1 AWT 

r4jUW ro!4W 


cluster 


.-.1,, C I \ 

cluster(ol) 


D07^ plAA 

B3 /o-G30y 


10818 


WrrLZ 


C/ICAilf CO 1 AWT 

F45UW rol4W 


cluster 


cluster (SI) 


B3oU-G313 










cluster(M) 


D/TJi pun 
B431-GUZ/ 










cluster (SI) 


D44J-1JUJ4 










cluster(SI) 


rS44fi-lJUj5 


10818 


WrrLZ 


C/1^A\X/ CQ1AWT 

r4jUW ro!4W 


cluster 


.-.1,, C FA 

cluster(M) 


B45 1 










possible star(I) 


d A ci niiJi 
B453-DU42 










cluster(^l) 


B4jfi-UU4y 


10407 


ALb/WrL 


F435W roUoW 


cluster 


cluster(^l) 


B4/5-V IZs 


10818 


\SICDA1 

WrrLZ 


C/1CAM7 CO 1 AWT 

F45UW rol4W 


cluster 


cluster (SI) 


B4oU-V ill 










cluster (SI) 


B483-D085 










cluster(SI) 


B484-G310 










cluster(SI) 


B486-G316 










cluster(S) 


B189D-G047 










cluster(SI) 


VDB0-B195D 


10818 


WFPC2 


F450W F814W 


cluster 


cluster(SI) 


NB21-AU5 


10006 


ACS/WFC 


F435W 


cluster 


cluster(SI) 


NB67 


10818 


WFPC2 


F450W F814W 


star 


star(SI) 


NB83 


5907 


WFPC2 


F555W F814W 


star 


star(SI) 
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Name 


Obs-ID 


Camera 


Filters 


Class HST 


Class C09 


B006D-D036 










cluster(SI) 


B012D-D039 










cluster(SI) 


B015D-D041 


10818 


WFPC2 


F450W F814W 


cluster 


cluster(SI) 


B111D-D065 


9794 


WFPC2 


F336W F439W F555W F675W F814W 


cluster 


cluster(SI) 


B206D-D048 










cluster(SI) 


B257D-D073 


10818 


WFPC2 


F450WF814W 


cluster 


cluster(I) 


DA047 










cluster(SI) 


V031 


10818(9709) 


WFPC2 


F450WF606W F814W 


cluster 


cluster(SI) 
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Table B.2. Classification of candidate young clusters listed in Tab. 2 of F05. 



Name 


Obs-ID 


Camera 


Filters 


Class HST 


Class C09 


rSUl J- VZU4 










cluster(M) 




6671 


WrrLZ 


rjjjW rol4W 


cluster 


cluster(Sl) 


B090 


10260 


ALi/WrL 


CAI\AWl UOl AW/ 


cluster 


cluster(SI) 


D1A1 i" 1 1 C A 

B101-G164 










cluster(SI) 


B102 


10260 


ACS/WrC 


F606W 


star 


star (SI) 


Bl 1 /-Gl /6 


9087 


WrrCz 


ci o £.w i 
r5 J6W 


cluster 


cluster(SI) 


B146 


101 18(5435) 


ACS/ WrC, WrrCZ 


1~7 1 ^ADHf CI C f CO AA\* T COl AWT 

rl60BW rzjjw rJOOW r814W 


cluster 


SLH 


B154-G208 


9087 


ACS/WFC 


F435W 


cluster 


cluster(SI) 


B 164-V253 










cluster(SI) 


Bly/-G247 










cluster(SI) 


B214-G265 










cluster(SI) 


B232-G286 


8059 


WhFLz 


OA/i\l/ Y2ACi\\\I C/' (\f Wl COl AW1 

riOOW F450W r606W r814W 


cluster 


cluster(SI) 


Dim pAi a 

Bzyz-GOio 


10631 


ACS/ WFC 


~c Al zwt TjeL(\cwi 

F435W r606W 


cluster 


candidate 


B311-G033 


6671(11081) 


WFPC2 


F555W F606W F814W 


cluster 


cluster(SI) 


B324-G051 


6699 


WFPC2 


F555W F814W 


cluster 


cluster(SI) 


B328-G054 


6699 


WFPC2 


F555W F814W 


cluster 


cluster(SI) 


B347-G154 


10818 


WFPC2 


F450W F814W 


cluster 


cluster(S) 


B423 


idate 








candidate 


B468 


5112 


WFPC2 


F555W F814W 


cluster 


cluster(I) 


NB16 


10818 


WFPC2 


F450W F814W 


cluster 


cluster(SI) 


B150D 










candidate 



